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Abstract 

Recent results are surveyed regarding the spectrum and eigenfunc- 
tions of the inverse square Calogero model with harmonic confinement 
and its relativistic analogue. 
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1 Introduction 

The Calogero model is a dynamical system that consists of particles on 
the line interacting pairwise through an inverse square potential and coupled 
to a harmonic external field. Both at the level of classical and quantum me- 
chanics the dynamics of the system has been studied in considerable detail 
in the literature Q. Key property of the Calogero model is that, despite 
the nontrivial interaction between the particles, the equations of motion de- 
scribing the dynamics can be solved in closed form. The exact solubility of 
the system stems from the fact that it is integrable, i.e., that there are as 
many independent integrals of motion (in involution) as degrees of freedom 
(viz. N). For the classical model the equations of motion were integrated by 
Olshanetsky and Perelomov using a Lax pair representation, whereas for the 
quantum system the spectrum and the structure of the eigenf unctions had 
already been determined before by Calogero 

More recently, Ruijsenaars and Schneider introduced a relativistic gen- 
eralization of the classical Calogero model without harmonic external field 
and solved the corresponding equations of motion 0, |§] . It was furthermore 
shown that also the relativistic system is integrable and that this integrability 
is preserved after quantization [^. At the quantum level, the Hamiltonian 
of the relativistic model is given by a difference operator rather than a dif- 
ferential operator; the nonrelativistic limit then corresponds to sending the 
step size of this difference operator to zero. 

Very recently the author has introduced a similar integrable relativistic 
analogue of the (quantum) Calogero model with harmonic external field [0], 
and computed the corresponding spectrum and eigenfunctions also for this 
case ip. The present contribution intends to provide an overview of these 
results as well as to describe some new developments that have led to a more 
explicit construction for the eigenfunctions of the nonrelativistic model ||^. 

2 The Quantum Calogero Model 

The quantum Calogero model is characterized by a Hamiltonian of the form 




(1) 
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It is immediate from the Hamiltonian that the model may be viewed as a sys- 
tem of N coupled harmonic oscillators. Should the inverse square coupling 
term be absent {g = 0), then it is of course very easy to solve the corre- 
sponding eigenvalue problem. The (boson) eigenfunctions are in that case 
the product of a Gaussian ground state wave function and symmetrized prod- 
ucts of Hermite polynomials, and the spectrum is that of an A^- dimensional 
harmonic oscillator. The remarkable observation by Calogero, to date around 
25 years ago, is that much of this picture is preserved after switching on the 
inverse square interaction {g > 0) [^f. 

Specifically, for the interacting system the ground state wave function 
becomes ^ 

^o(^) = n \xj - Xkf exp{-- ^i) (2) 

i<j<k<N ^ l<j<N 

and the ground state energy reads 

Eo = ujN{l + {N-l)g). (3) 

The wave functions of the excited states are again products of the ground 
state wave function and certain symmetric polynomials: 

^^(f) = ^o(x) Pnix), (4) 

where n = (ni, . . . , n^) denotes a vector of (integer) quantum numbers with 
ni > n2 > ■ ■ ■ > riN > labeling the eigenfunctions. The corresponding 
eigenvalues are given by 

Eh = Eo + 2uj ^j- (5) 

It is clear from Eq. ^ that, apart from an overall shift in the energy, the 
spectrum of the model with inverse square interaction coincides with that 
of a system of independent harmonic oscillators. Crucial difference with the 
latter system, however, is that for the interacting system the wave functions 
"^nix) do not separate in one-particle wave functions. In fact, until recently 
not much information had been available regarding the precise nature of the 
polynomials Pnix), except for small particle number. (For < 5 the poly- 
nomials Pfi{x) were constructed explicitly by Perelomov and Gambardella 

MM-) 
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3 Creation and Annihilation Operators 



It is well known that the wave functions for a system of independent har- 
monic oscillators can be constructed by means of creation and annihilation 
operators. Interestingly enough, it was recently discovered that the same 
classical technique may also be applied to the case with an inverse square 
coupling between the oscillators 

To this end it is convenient to introduce the concept of the so-called Dunkl 
derivative | |12| ]: 

d 

= ^ g Yj (^k- xjy^Sj^k, (6) 



dXi 



where Sj^k denotes the transposition operator interchanging the particles j, 
k 



Sj- ...,Xj,...,Xk,..., Xn) = '^{Xl, ...,Xk,...,Xj,..., xn)- (7) 

The Dunkl derivative Dj is a deformation of the ordinary derivative dj = 
d/dxj with the coupling constant g acting as a deformation parameter. The 
essential point is now that it is possible to incorporate the inverse square 
coupling between the oscillators in the standard construction of the eigen- 
functions by means of creation and annihilation operators if one replaces the 
ordinary derivative by the Dunkl derivative. More specifically, after setting 



-^{tDj + uXj] 



(8) 



one has 



Hc^ E (A^A- + A-A^) 



and 



Hc,Af 



±2ujA^. 



+ E 



dig- Sj, 



Xkj 



(9) 



(10) 

When restricted to the boson sector (i.e., the space of permutation invariant 
wave functions), the operator He coincides with the Calogero Hamiltonian 
He (cf. Eq. IID. The ground state wave function in this sector is deter- 
mined by its permutation-invariance and the fact that it is annihilated by 
the lowering operators 



AJ ^oix) 



0, 



J 



N. 



(11) 
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The excited states are obtained by acting with the creation operators on the 
ground state 

^n^(^) = E ( n (^^0))"0*0(X) (12) 

aes^ i<j<N 

(where we have, in order to obtain boson eigenfunctions, symmetrized over all 
permutations). The shift of the spectrum as compared to the model without 
inverse square coupling between oscillators originates from a small change 
in the usual commutation relations satisfied by the modified creation and 
annihilation operators A^: 

[Af,A^] = 0, [Aj,At]=ujil + g ^ S,,i)6,,k - ujgS,,,. (13) 

l</<Ar 

(So i/c^o = Ej{[Aj,Aj] + 2A+A7)^o = Ejuj{l + {N- l)g)^o = ^o^o, 
with Eq given by Eq. ||.) 

4 Relativistic Analogue 

A few years ago, Ruijsenaars introduced an integrable quantum iV-particle 
system characterized by a (rather unorthodox) Hamiltonian given by the 
second order difference operator 

Hr= (^/^' + e'"^ yl^^) , (14) 

where exp(±r^)^(xi, . . . , xn) = *(a;i, . . . , Xj-i,Xj =F i, Xj+i, . . . xn) and 

Vj = n v{xj-Xk), v{z) = 1 + g/{iz) (15) 

l<fe<iV, k=ij 

(with Vj denoting the complex conjugate of Vj). As it turns out, Ruijsenaars' 
difference model may be interpreted as a system composed of N relativistic 
particles in [1 + 1) D that interact with each other by means of the coefficients 
Vj. For g = the particles are independent (Vj = 1), whereas for g > 
each particle feels the presence of the remaining — 1 particles as a change 
of its dynamical mass. A more detailed discussion of the model (both at 
the classical and quantum level), with an emphasis on matters involving 
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integrability and Poincare invariance, can be found in Ruijsenaars' papers 
i 1,1- 

Quite recently, the present author observed that it is possible to introduce 
an external field coupling to the relativistic system without destroying its 
integrability 0. The Hamiltonian of the corresponding difference model 
reads 

H= i^l'^ ' + ' e"^ r//' - Vj - , (16) 

with 

Vj = w{xj) Y[ v{xj-Xk), (17) 
v{z) = l + g/{iz), w{z) = {a + iz){b + iz). (18) 

The function w encodes the external field. For w = 1 (this is achieved by 
sending the parameters a and b to infinity after having rescaled the Hamilto- 
nian by division by ah), the Hamiltonian H reduces, up to an irrelevant ad- 
ditive constant, to the Ruijsenaars Hamiltonian H^i. The constant is caused 
by the part —J2i<j<N(yj + Vj), which does not depend on Xj ii w = 1 (as 
is readily seen with the aid of Liouville's theorem after having inferred that 
the expression is regular in xj and bounded for xj oo). 

Below, we will discuss the spectrum and eigenf unctions of the difference 
model with external field given by the Hamiltonian in Eqs. p^|T8| and describe 
its relation to the Calogero model discussed in Sections |^ and |[ Throughout 
it will be assumed that the coupling constant g be nonnegative and that (the 
real parts of) the parameters a and b be positive. 



5 Spectrum and Eigenfunctions 

The ground state wave function for the difference Hamiltonian H is given by 



n 

i<j<k<N 



T{g + i{xj -Xk)) 



T{i{xj - Xk)) 



Y[ \T{a + iXj)T{b + iXj)\. (19) 



Probably the simplest way to see that this is indeed an eigenfunction is to 
check that conjugation of the Hamiltonian with \E'o yields (using the standard 
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difference equation r{z + 1) = zr{z) for the gamma function) 



a 



n = ^o^^^o = i^J i^"'^' - 1) + (e - 1)) • (20) 

Tlie transformed operator Ti. clearly annihilates constant functions, so \E'o is 
an eigenfunction of H with eigenvalue zero. 

Just as for the nonrelativistic Calogero model, the wave functions corre- 
sponding to the excited states are a product of the ground state wave function 
and symmetric polynomials. More precisely, one has 



with 



and 



H^n = Eh (21) 



En= Y Ujirij + a + a + b + b + 2{N - j)g) (22) 



^h{x) = ^o{x)Ph{x), (23) 

where Pfi{x) is the symmetric polynomial determined by the conditions: 

i. Ph{x) = nifiix) + Y ^n,n' rufi'ix); 

n'<n 

/oo roc 
■ ■ ■ / Pnix) mfi'{x) "^lix) dxi ■ ■ ■ dxN = if n' < n. 
-oo J —oo 

Here the functions mf{{x), n = (ni, . . . , un) with rii > n2 > ■ ■ ■ > un > 0, 
denote the basis consisting of symmetrized monomials 

mn{x) = J2 Xi^-'-x''^'', (24) 

n'e5]v(n) 

which is partially ordered by the definition 

n' <n iff n[-\ h n'^ < ni H h for k = l,...,N (25) 

{n' < n \S. n' < n and n' ^ n). Thus, the polynomial Pn{x) boils down 
to the symmetrized monomial m^(x) minus its orthogonal projection with 
respect to the inner product with weight function ^'o(^) onto the finite- 
dimensional subspace spanned by the monomials mfi'{x) with n' < n. For 
= 1 the resulting polynomials are well-studied in the literature and known 
as continuous Hahn polynomials |T3|, P|. 
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The proof of the above statements hinges on a standard technique going 
back (essentially) to Sutherland |15]. First it is shown that the transformed 
operator Ti is triangular with respect to the partially ordered monomial basis, 
i.e. 

{Hmji){x) = J2 [^]n,n' nT'H'ix), (26) 

n'<n 

where the [7i]^,,j' represent certain (complex) matrix elements. It is clear that 
acting with 7i on a monomial yields a permutation invariant rational 
function. The permutation symmetry guarantees that the simple poles at 
Xj = Xk (caused by the zero in the denominator oiv{z)) all cancel each other. 
Hence, the rational function is actually a (permutation invariant) polynomial 
and can thus indeed be expanded in symmetrized monomials. To see that 
in this expansion only monomials m^/ with n' <n occur (triangularity), one 
uses the asymptotics at infinity. Setting Xj = K^^ with yi > y2 > ■ ■ ■ > Un > 
gives the following asymptotics for R +00: 

m^/ = R^'-y +o{R^'y), (27) 

HrriH = 0{R^-y). (28) 

By comparing the asymptotics of Eqs. ^ and and using the fact that 

n' <n iff n' ■ y < n - y with yi > y2 > ■ ■ ■ > yN > 0, (29) 

one infers that the matrix elements [?^]n,n' in the expansion of Tirufi in terms 
of symmetrized monomials m^' can only be nonzero if n' < n, which gives 

Next one observes that since H is Hermitian the transformed operator 
7i is symmetric with respect to the inner product with weight function 
^lix): 

/OO roo 
• ■ ■ / {Hm^){x) mff{x) ^o(^) '^^1 ■ ■ ■ dxN = 
-00 J —00 

/OO rco 
• • • / m^(x) {Hmfii){x) "^lix) dxi ■ ■ - cLxn. (30) 
-00 J —00 

By combining the triangularity and symmetry of Ti. it is not difficult to 
deduce that T-CPa is a linear combination of monomials m^/ with n' < n, 
which is orthogonal to m^' with n' < n. Hence, TiPft must be proportional 
to Pfi, i.e., Pfi is an eigenfunction of H. The corresponding eigenvalue Efi is 
obtained by computing the diagonal matrix element [7i]ra,n, i-e., by explicitly 
computing the leading coefficient in Expansion E6|. 
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6 The Nonrelativistic Limit 

To relate the difference model to the Calogero model one needs to explic- 
itly introduce a step size parameter by rescaling the positions xj and the 
parameters a, b. Substituting Xj P~^Xj (so dj Pdj), a {P'^uj)~^, 
b — >■ (/3^a;')~^, and multiplying H by (3'^ujuj' , leads to a Hamiltonian H given 
by Eqs. |16|, |T7|with exp(±ic}j) replaced by exp(±i/99j) and functions f , w of 
the form 

v{z) = 1 + (3g/{iz), w{z) = (3~'^{l + i(3ujz){l + i(3uj' z). (31) 

The step size parameter (5 should be compared with the inverse of the light 
speed appearing in Ruijsenaars' model 0. 
For /5 — > one now has 

«^ Y.{-^. + f-f) + E 77^-^0. (32) 

with £o = {oj + u')N{l + {N —l)g), dj^d Efi ^ 2{ijj + uj')Y.i<j<N''^j- The wave 
functions go (after dividing by a divergent numerical factor arising from the 
gamma factors in \E^o; Eq. |l^) over in ^'^(x) = \E'o(x)P^(x), where 

^o(x)= n \^j-^k\' n e-^(-+-')^? (33) 

^<j<k<N l<j<N 

and Pn{x) is the polynomial determined by the Conditions i. and ii. of 
Section ^ with now taken from Eq. This way we recover for (3^0 
the Hamiltonian, the spectrum and the eigenfunctions of the Calogero model. 
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